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We present first-principles calculations on electron transport through Na nanowires at finite bias 
voltages. The nanowire exhibits a nonlinear current-voltage characteristic and negative diflterential 
conductance. The latter is explained by the drastic suppression of the transmission peaks which is 
attributed to the electron transportability of the negatively biased plinth attached to the end of the 
nanowire. In addition, the finding that a voltage drop preferentially occurs on the negatively biased 
side of the nanowire is discussed in relation to the electronic structure and conduction. 
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Currently, quantized conductances of atomic 
nanowires with a few atoms width and nanometer 
length, which are in steps of Go = 2e^//i (e: the 
electron charge, h: Planck's constant), have attracted 
much attention because of both fundamental and 
technological importance. A number of exciting ex- 
perimental and theoretical studies on the conductance 
quantization for various metal nanowires have been 
reported§ | | |, |, |, | g 0, P. Above 
all, the studies for Na nanowires arc of particular 
interest, because the nature of electron conduction 
through nanowires is expected to emerge best for Na 
nanowires having the simplest electronic structure 
where monovalent electrons predominate. Elaborate 
investigations on electron-transport properties of Na 
nanowires are, however, not many yet. A few exper- 
iments have been reported^, where the variation of 
the conductance was measured during the withdrawal 
of the Na film to show the stairlike behavior of the 
quantized conductance. Theoretically, there have been 
some first-principles calculations fully self-consistently 
performed using the model incorporating a couple of 
semi-infinite electrodes so as to evaluate electron flow 
through Na nanowires ||, ^, 0]- But, these theoretical 
studies are limited to zero bias voltage |^, ^ or the low 
bias regime 0. The analysis and discussion of electron 
transport through the Na nanowire should be extended 
to a higher bias voltage in order to gain insight into 
a unique conduction property possible under a strong 
electric field. 

In this letter, we present a theoretical analysis of elec- 
tron transport through the Na nanowire under the ap- 
plication of a sizable bias voltage, based on the den- 
sity functional theory within the local density approxi- 
mation. All the calculations determining the electronic 
structure and transport properties are self-consistently 
evaluated by the method making use of the Lippmann- 
Schwinger equation proposed by Lang|^. The norm- 
conserving pseudopotentials||l^ are adopted to describe 
the electron- ion interaction, and the real-space timesav- 
ing double-grid technique |l^ is used by virtue of an effi- 
cient computation of the nonlocal parts of the pseudopo- 
tentials. We will demonstrate that the Na nanowire ex- 



hibits an unusual nonlinear current-voltage (I-V) charac- 
teristic which leads to negative differential conductance 
behavior. Notice that the negative differential conduc- 
tance is expected to be an essential property to real- 
ize fast switching in future electronic devices. This pa- 
per is probably the first report describing that a sim- 
ple nanowire in which all atoms are tightly binding to- 
gether possesses the negative feature of the differential 
conductance [0. We will also discuss how an applied 
bias voltage gives rise to the change of the effective poten- 
tial along the nanowire (voltage drop) and of the charge 
density distribution. 

According to Refs. || and 0, we here adopt the model 
composed of three parts for the nanowire system: a 
straight atomic nanowire, two plinths, and a couple of 
semi-infinite jellium electrodes sandwiching them (see 
Fig. gfor illustration). The straight nanowire is made up 
of N (Af=l~5) Na atoms equispaced at 7.0 a.u., which 
is equal to the nearest neighbor atomic distance of the 
Na crystal. Each plinth is modeled after the Na(OOl) 
ideal surface, i.e., it is formed by four Na atoms, one at 
each corner of a square with the side length of 8.1 a.u. 
These plinths are attached to the two ends of the straight 
nanowire and act as buffer layers to moderate the un- 
physical effect due to an oversimplified electrode model 
of jellium. The electron density of the jellium electrodes 
is taken equal to the average valence-electron density of 
the Na crystal (r^ = 4.0 a.u.)j^. The distance between 
the jellium edge and the atoms forming the plinth is set to 
be 4.0 a.u. In the present calculation, we do not consider 
any structural relaxation. 

In the numerical calculations, we employ the super- 
cell on which periodic boundary conditions are imposed 
in the directions perpendicular to the nanowire {x and y 
directions) and nonperiodic boundary condition in the di- 
rection parallel to it (z direction). Therefore, wave func- 
tions are expressed by the Laue representation. The side 
length of the supercell in the x and y directions is 24.0 
a.u. and the number of the plane waves in these direc- 
tions is taken to be 24^, which corresponds to a cutoff 
energy of 7.4 Ry. We use only the T point in the xy 
Brillouin zone because of a large supercell adopted here. 
The grid size in the z direction is chosen to be 0.5 a.u. 
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FIG. 1: (a) Conductance and (b) LDOS around the straight 
nanowire of the respective systems at the zero bias limit. In 
(b), the broken line represents the Fermi level. The LDOS 
curves are displaced vertically for clarity. 



Before moving on to the discussion of electron- 
transport properties for the Na nanowire system at a 
finite bias voltage, we briefly report the results obtained 
at the zero bias limit. In Fig. |l](a), the calculated con- 
ductance is drawn as a function of the number of atoms 
N in the straight nanowire: in the case of iV = 1, i.e., 
the case of the closed packed structure of Na atoms, the 
conductance takes the value 2.860, while in the cases of 
N =2~5, it is close to the quantized value of IGqI^. 
These features of the conductance are in good agreement 
with the experimental results obtained by Krans et at, 
where the conductance histogram shows the three dis- 
tinct points: a peak at slightly lower than the 3Go quan- 
tized value, a peak at IGq and a missing 2Go peak@. 
In Fig. ||(a), an oscillatory behavior in the conductance 
is observed for A'^ = 2 5, i.e., the even-atom-number 
nanowire systems exhibit slightly lower conductance than 
the quantized value of IGq, and the odd-atom- number 
systems show fully quantized conductance IGq. This os- 
cillation can be understood in terms of the local density 
of states (LDOS) around the nanowire at the Fermi level: 
as seen in Fig. |l|(b), the LDOS for the even-atom- number 
(odd-atom-number) nanowire has an off-resonance (on- 
resonance) transport property at the Fermi level [|l5| . 
Based on the Friedel sum rule, Sim et al. recently pre- 
dicted the occurrence of the even-odd oscillation for the 
Na nanowire sandwiched by the pyramidal-shape tips 
I prf . Our present results verify their prediction of the 
even-odd oscillation. 

We now discuss how electron-transport properties of 
the nanowire are affected when a sizable bias voltage is 
applied. The range of bias voltages is 0.0~3.5V and a 
long nanowire of iV = 5 is adopted. The current per 
supercell is simply expressed as the integration of the 
transmission over the voltage window, which represents 
the energy range between the Fermi levels of the two elec- 
trodes. Hereafter, the respective electrodes are called left 
and right electrodes, and the bias voltage is applied with 
the left electrode negative and the Fermi level of the right 
electrode kept fixed. Figure |2| shows the I-V characteris- 
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FIG. 2: I-V characteristic of the N = 5 Na nanowire system 
(solid line) and its differential conductance (broken line). 
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FIG. 3: (a) Transmission and (b) LDOS around the straight 
nanowire of electrons incident from the left electrode. The 
energy is measured from the Fermi level of the right electrode. 
The shaded areas show the voltage windows. The eight curves 
in each panel correspond to the bias voltages indicated, and 
are displaced vertically for clarity. 



tic of the system and its differential conductance. Here, 
the differential conductance is defined as the derivative of 
the current with respect to the bias voltage. It is clearly 
seen that the system exhibits the nonlinear 1-V character- 
istic and in particular, the differential conductance shows 
unusual behavior near zero bias voltage and around 2.5 
V; the differential conductance rapidly decreases in the 
range of 0.0~0.5 V, deviating from the quantized value 
IGo observed at the zero bias limit, and moreover, it is 
negative at the bias voltages of 2.3~2.7 V with a mini- 
mum at 2.5 V. 

We first examine why the reduction of the differential 
conductance is induced at the low bias voltages of 0.0'--^0.5 
V. The transmission and the LDOS around the straight 
nanowire of electrons incident from the left electrode are 
depicted as a function of the electron energy measured 
from the Fermi level of the right electrode in Figs. ||(a) 
and (b), respectively. The right edges of the voltage win- 
dows (shaded areas) in these figures correspond to the 
Fermi level of the left electrode for the respective bias 
voltages. When the bias voltage is boosted from 0.0 V to 
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FIG. 4: LDOS around the left plinth at the bias voltages 
of 0.0, 2.0, 2.5 and 3.0 V. The energy is measured from the 
Fermi level of the right electrode. The curves are displaced 
vertically for clarity. 



0.5 V, a dip in the transmission curve, which is indicated 
by D in Fig. ^(a), appears in the voltage window, caus- 
ing the reduction in the differential conductance. This 
transmission dip D is brought about by the off-peak of 
the LDOS marked D in Fig. |(b) 

Next, we discuss the negative feature of the differential 
conductance at around 2.5 V seen in Fig. ^. In the trans- 
mission in Fig. ^(a), one can see that the leftmost peak 
clearly observed at the bias voltage of 2.0 V is suddenly 
suppressed at 2.5 V (A^A'), and that one more peak at 
2.5 V is also suppressed when the bias voltage is raised to 
3.0 V (B^B')||l^. As a consequence of the suppression, 
the integration of the transmission over the voltage win- 
dow decreases upon raising the bias voltage from 2.0 V 
to 3.0 V, which leads to the negative differential conduc- 
tance. In contrast, such a disappearance of the transmis- 
sion peak is no longer observed at the bias voltage of 3.5 
V. It should be noticed that the suppression of the trans- 
mission peaks is in accordance with the vanishing of the 
corresponding peaks in the LDOS around the straight 
nanowire of electrons incident from the left electrode, 
marked by A-^A' and B-^B' in Fig. |^(b). This van- 
ishing implies that some electrons coming from the left 
electrode with specific energies cannot reach the nanowire 
when the bias voltage exceeds a certain threshold. 

In order to investigate why these transmission peaks 
and LDOS peaks around the straight nanowire suddenly 
disappear as the bias voltage is raised, we focus on the 
electronic structure around the left plinth. Figure ^ shows 
the LDOS around the left plinth of electrons incident 
from the left electrode at the bias voltages of 0.0, 2.0, 
2.5 and 3.0 V, where each curve shifts to the higher en- 
ergy side according to the applied bias voltage. This shift 
corresponds to the linear rising of the effective potential 
around the left plinth upon applying bias voltage, which 
will be discussed later for Fig. |^. We see in Fig. ^ that 
the leftmost LDOS peaks appear at 0.0, 0.5 and 1.0 eV 
at the bias voltages of 2.0, 2.5 and 3.0 V, respectively. 
It is now common knowledge that the energy of the left- 
most LDOS peak corresponds to the first onset energy 
of the transmission 1^, |^, ||. Consequently, when the bias 
voltage of 2.0 V is applied, electrons with energy lower 
than 0.0 eV are refused to go through the left plinth. 
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FIG. 5: (color) Shifts of the local effective potential along 
the axis of the nanowire upon applying the bias voltages of 
2.0, 2.5, 3.0 and 3.5 V. The atomic geometry is shown in the 
panel. The Na atoms forming the N = 5 straight nanowire 
are labeled 1~5. 
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FIG. 6: (color) Changes in the charge density distribution 
upon applying the bias voltages of (a) 2.5 V and (b) 3.5 V. 
The red shows the increase and the blue the decrease. An 
area with sides of 55.0 a.u. x 23.0 a.u. is displayed. The Na 
atoms forming the N = 5 straight nanowire are labeled 1~5. 
Broken lines represent the edges of the jellium electrodes. 



while those with energy above 0.0 eV readily enter the 
plinth and the nanowire to yield the transmission peak 
A in Fig. ^(a). In the case of the bias voltage of 2.5 V, 
the threshold energy for electrons to pass through the left 
plinth is 0.5 eV, therefore, electrons with the energy of 
0.0 eV are prohibited from entering the left plinth. Thus, 
as the bias voltage is boosted from 2.0 V to 2.5 V, the 
suppression of the transmission peak at 0.0 eV, A^A', 
is brought about. This mechanism of the suppression of 
the transmission peak also holds in the case of B^B' in 
Fig.|(a). 

Finally, we explore how the applied bias voltages drop 
along the nanowire. Figure || shows the differences of the 
local effective potential along the axis of the nanowire 
at some significant bias voltages from that at zero bias 
voltage, VeffiVB) - K//(0.0), at Vb =2.0, 2.5, 3.0 and 
3.5 V. In Fig. @, one can see that the intensive volt- 



4 



age drop takes place mainly around the negatively bi- 
ased end of the nanowire (atoms numbered 1 and 2 in 
Fig. while the voltage drop is nominal near the 

two plinths and on the positively biased side (atoms 
3~5). The localization of the voltage drop is understood 
as follows. Electrons abundant in proportion to the ap- 
plied bias voltage rush into the left end of the nanowire 
from the left electrode through the plinth. The conges- 
tion of electron flow occurs around the region between 
atoms 1 and 2, where the narrowing of the path along 
which electrons pass takes place due to the geometri- 
cal transition from the wider electrode to the narrower 
nanowire. Although only electrons with some specific 
momenta can enter the conduction channel and reach 
the opposite electrode without any backscattering, the 
remaining electrons suffer from a considerable amount of 
reflection. Thus, the resistance for electron flow localizes 
around the region between atoms 1 and 2, which results 
in the intensive voltage drop there. This mechanism of 
electron flow is analogous to the vena contracta effect 
known in hydrodynamics. The absence of a sizable volt- 
age drop near the two plinths and on the positively biased 
side of the nanowire is associated with the change in the 
charge density distribution upon applying bias voltage. 
Figure ^ shows the profiles of the difference of the charge 
densities, p{Vb) - p(O.O), at Vb =2.5 and 3.5 V. On the 
right of the atoms forming the left plinth and the left of 
atom 2, one observes considerable increases of the charge 
density. These increases reflect the screening effect which 
keeps the electric potential of the negatively biased jel- 



lium electrode and plinth equal to the negative electrode 
potential applied to the left jellium. On the other hand, 
the most significant decreases of electrons occur on the 
right of atom 2 so that the electric potential of the pos- 
itively biased side of the nanowire system (atoms 3'^5, 
the right plinth and jellium electrode) becomes close to 
the positive electrode potential assigned to the right jel- 
lium. As a consequence of the negative electric potential 
on the left of atom 1 and the positive one on the right 
of atom 2, the localization of the voltage drop between 
atoms 1 and 2 is naturally realized. 

In conclusion, a first-principles calculation method 
using the Lippmann-Schwinger equation was imple- 
mented to elucidate electron-transport properties of Na 
nanowires at finite bias voltages. Employing a nanowire 
model including the plinths as well as the semi-infinite 
jellium electrodes, we obtained the nonlinear I-V charac- 
teristic and the unusual behavior of the differential con- 
ductance. We demonstrated that the negative feature of 
the differential conductance is explained by the suppres- 
sion of the transmission peaks due to the electron trans- 
portability of the negatively biased plinth. It was also 
found that the voltage drop preferentially occurs around 
the negatively biased end of the nanowire, and this find- 
ing was discussed in relation to the electronic structure 
and conduction of the nanowire. 

This work was supported by a Grand-in-Aid for COE 
Reseach (No. 08CE2004) from the Ministry of Education, 
Culture, Sports, Science and Technology. 



[1] J. M. van Ruitenbeek in Metal Clusters at Surfaces Struc- 
ture, Quantum Properties, Physical Chemistry, edited by 
K. H. Meiwes-Broer (Springer, Berlin, 2000), and refer- 
ences therein. 

[2] J. M. Krans, J. M. van Ruitenbeek, V. V. Fisun, I. K. 
Yanson, and L. J. de Jongh, Nature (London) 375, 767 
(1995); A. I. Yanson, I. K. Yanson, and J. M. van Ruiten- 
beek, Nature (London) 400, 144 (1999). 

[3] H. Ohnishi, Y. Kondo, and K. Takayanagi, Nature (Lon- 
don) 395, 780 (1998). 

[4] N. D. Lang, Phys. Rev. Lett. 79, 1357 (1997); Phys. Rev. 
B 55, 4113 (1997). 

[5] A. Nakamura, M. Brandbyge, L. B. Hansen, and K. W. 
Jacobsen, Phys. Rev. Lett. 82, 1538 (1999). 

[6] N. Kobayashi, M. Brandbyge, and M. Tsukada, Phys. 
Rev. B 62, 8430 (2000); Surf. Sci. 433-435, 854 (1999). 

[7] S. Tsukamoto, Y. Fujimoto, T. Ono, K. Inagaki, H. Goto, 
and K. Hirose, Mater. Trans. 42, 2253 (2001). 

[8] J. C. Cuevas, A. LevyYeyati, A. Martm-Rodero, G. R. 
Bollinger, C. Untiedt, and N. Agrait, Phys. Rev. Lett. 
81, 2990 (1998). 

[9] N. D. Lang, Phys. Rev. B 52, 5335 (1995). 
[10] R. N. Barnett and U. Landman, Nature (London) 387, 
788 (1997); M. Okamoto and K. Takayanagi, Phys. Rev. 
B 60, 7808 (1999). 



[11] H.-S. Sim, H.-W. Lee, and K. J. Chang, Phys. Rev. Lett. 
87, 096803 (2001). 

[12] N. Troullier and J. L. Martins, Phys. Rev. B 43, 1993 
(1991); K. Kobayashi, Comput. Mater. Sci. 14, 72 (1999). 

[13] T. Ono and K. Hirose, Phys. Rev. Lett. 82, 5016 (1999). 

[14] The negative differential conductance was found for sys- 
tems in which a pair of atomic clusters, each coupled to 
an electrode, are weakly interacting at a distance from 
each other, e.g., by N. D. Lang, Phys. Rev. B 55, 9364 
(1997), and by I.-W. Lyo and Ph. Avouris, Science 245, 
1369 (1989). 

[15] N. D. Lang and P. Avouris reported a similar oscillatory 
behavior and its interpretation for carbon nanowires in 
Phys. Rev. Lett. 81, 3515 (1998). 

[16] The spacial distribution of the LDOS relevant to the peak 
A (B) shows an eigenmode pattern with a loop (node) 
at the center of the nanowire. Because of this distinct 
difference, the peak B is not regarded as a shifted one of 
the peak A to the higher energy side. 

[17] A similar behavior of the voltage drop was obtained by 
M. Brandbyge, N. Kobayashi, and M. Tsukada using a 
parameterized tight-binding method for Au nanowires, 
where s- valence electrons also predominate as in the case 
of Na nanowires. See Phys. Rev. B 60, 17064 (1999). 



